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Abstract 

Our previous calculations of the sea- and valence-quark mass dependence of the 
pseudoscalar meson masses and decay constants is repeated on a 16^ • 32 lattice which 
allows for a better determination of the quantities in question. The conclusions are 
similar as before on the 16^ lattice pp. The two light dynamical quark flavours we sim- 
ulate have masses in the range jrug < mu,d < ^rug. The sea quark mass dependence 
of /tt and m^/mq is well described by the next-to-leading order (NLO) Chiral Pertur- 
bation Theory (ChPT) formulas and clearly shows the presence of chiral logarithms. 
The valence quark mass dependence requires the presence of NNLO contributions in 
Partially Quenched ChPT (PQChPT) - in addition to the NLO terms. The 0{a) 
lattice artifacts in these quantities turn out to be small. 
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1 Introduction 



In Quantum Chromodynamics (QCD) - the theory of strong interactions - there are 
two very hght quarks and one moderately Hght quark (u, d and s, respectively). The 
strong interaction dynamics at low energies can be formulated by an effective chiral 
Lagrangian which incorporates the symmetry constraints following from the sponta- 
neously broken chiral symmetry of the light quarks. In this low energy effective theory 
the interactions are described by a simultaneous expansion in powers of momenta and 
light quark masses [21 El • The coefficients of the interaction terms in the effective chi- 
ral Lagrangian - the Gasser-Leutwyler constants ~ are free parameters which can be 
constrained by experimental data and also calculated from the underlying basic QCD 
Lagrangian in the framework of the non-perturbative lattice regularization. 

In numerical lattice QCD simulations the quark masses are free parameters. Chang- 
ing these parameters gives an excellent opportunity to precisely determine the Gasser- 
Leutwyler constants. In fact, ChPT based on the chiral Lagrangian can be extended 
by changing the valence quark masses in quark propagators independently from the 
sea quark masses in virtual quark loops. This leads to partially quenched chiral per- 
turbation theory (PQChPT) 01. 

The aim of numerical simulations in QCD is to reach the regime of light quark 
masses where NLO chiral perturbation theory gives a good approximation. In previ- 
ous papers |Slinm our collaboration started a series of simulations with two equal-mass 
light quarks {qq) with the goal of extracting the values of the Gasser-Leutwyler con- 
stants conventionally denoted by L^, (k = 1,2,...). Later on it will be possible to 
extend these calculations by also including the s-quark (qq+q). 

In our previous paper we started some larger scale simulations on a 16^ lattice at 
the gauge coupling /3 = 5.1 which corresponds to a lattice spacing of a ~ 0.2 fm. Since 
it became clear that interesting results can be obtained already at this relatively rough 
discretization scale, we decided to repeat and extend these simulations on a IG'^ • 32 
lattice which is better suited for extracting quantities like the pseudoscalar ("pion") 
mass (m,r) and decay constant (/-n-). Our work profited from the valuable experience 
of previous simulations by other collaborations [71|Slini • 

Since the present work is on the same topics as ref. 0, we shall often only refer to 
it without repeating its full content. In general, we use the conventions and notations 
of ^EllEl- Nevertheless, we also try to make the present paper easily understandable 
for the reader and therefore repeat the main definitions and relations. In the next 
section we deal with the sea quark mass dependence of and ttZtt. In section O the 
valence quark mass dependence is considered and the question of the magnitude of 
leading lattice artifacts is investigated. Section jl] is a short summary of our experience 
with the Monte Carlo updating algorithm. The last section contains the summary and 
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discussion. 



2 Sea quark mass dependence 

We performed Monte Carlo simulations with Ng = 2 degenerate sea quarks on a 163-32 
lattice at gauge coupling /3 = 5.1 and four values of the hopping parameter k: kq = 
0.176, Ki = 0.1765, K2 = 0.1768 and K3 = 0.177. Three of these points have also been 
simulated previously on the 16^ lattice in ref. pp. The point at ki = 0.1768 is new. We 
collected 950-1000 gauge configurations per point which are typically separated by 10 
update cycles consisting out of boson field and gauge field updates and noisy correction 
steps. (Some observations about the algorithm will be summarized in section |2) 

A collection of the values of some basic quantities in these simulation points is 
given in tabled the Sommer scale-parameter in lattice units ro/a, the pion mass in 
lattice units auiT^, the quark mass parameter Mr = (ro?n.7r)^, the bare PCAC quark 
mass Zqaniq including the multiplicative renormalization factor Zq = Zp/Za, the ratio 
of the PCAC quark masses with respect to the reference sea quark mass at k = 
and the pion decay constant in lattice units o/tt divided by the renormalization factor 
Z^. (The normalization of the pion decay constant is such that the physical value is 
U ^ 93MeV.) 



Table 1: The values of some basic quantities in our simulation 
points. Statistical errors in last digits are given in parentheses. 









1^2 




ro/a 


2.229(63) 


2.212(44) 


2.621(46) 


2.528(51) 




0.6542(10) 


0.5793(17) 


0.3919(46) 


0.3657(24) 


Mr = {rom^y 


2.13(12) 


1.642(72) 


1.055(36) 


0.855(34) 


Zqaniq 


0.07092(27) 


0.05571(30) 


0.02566(27) 


0.02208(21) 


CTi = rUqi/rnqQ 


1.0 


0.7856(56) 


0.3618(44) 


0.3113(31) 




0.2819(15) 


0.2590(14) 


0.2008(17) 


0.1936(16) 



Comparing tabled to the corresponding one (table 3) in 01 one can see that these 
quantities extracted on the 16^ • 32 lattice differ considerably from those extracted on 
the 16^ lattice. The change of ro/a is about 2-5%. The difference in am^ increases 
from 3% at kq to about 16% at K3 whereas Zqamq differs at kq by 5% and at K3 
already by about 28%. However, as we shall see later on, considering ratios of the pion 
mass-square and of the pion decay constant as a function of the ratios of PCAC quark 
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masses (denoted by a for sea quark masses and ^ for valence quark masses) it turns 
out that almost all changes between the 16^ and 16^ • 32 lattices cancel. 

In tabled the bare quark mass obtained from the PCAC relation is shown: niq = 
^PCAC ^ (For details of its numerical determination see section 3.1.1 in ref. [Hj.) An- 
other possibility to define the quark mass is to take arriren = tJ-ren = ■^m(/Uo — fJ-cr) 
where /io = 1/(2k) — 4 is the bare quark mass in the Wilson- fermion action, /Ucr is its 
critical value corresponding to zero quark mass and is an appropriate multiplica- 
tive renormalization factor. The values of corresponding to k,q,...,K3 are /xo{o) = 
-1. 1590909.., /io{i) = -1.1671388.., //o(2) = -1.1719457.., /xo{3) = -1.1751412.., re- 
spectively. Comparing the values of Zqaruq or ai in tabled to the values ^o(j) o^^^ can 
see that the relation between them is highly non-linear. This implies the same also for 
the relation between the (ratios of) rriq^^^ and rriren- The non-linear terms in this 
relation are lattice artifacts which have to vanish in the continuum limit but they are 
large at our lattice spacings. 

A consequence of the strongly non-linear relation between a and /io is that the 
determination of /Xcr (or Hcr) has large uncertainty. In fact, with our four points 
only we could not find a convincing extrapolation of a to zero. A crude quadratic 
extrapolation gives iicr = —1.180(4) or Kcr = 0.1773(2). The uncertainty in the critical 
point implies an uncertainty in the extrapolation of physical quantities, too, which is 
necessary in a quark mass independent renormalization scheme. In case of the lattice 
spacing, which can be obtained from the extrapolation of tq/g to the critical point, 
table n shows that the values of rg/a increase between kq and K2 but between K2 
and K3 they are within errors constant. Therefore we take this constant value as the 
extrapolated one: [ro/a]cr = 2.57(5). This gives, with tq = 0.5 fm, for the quark mass 
independent lattice spacing a = 0.195(4) fm. 

In the ChPT formulas the quark mass can be represented by the dimensionless 
quantity 

_ 2Bomg 

JO 

where Bq is a conventional parameter with dimension mass and /q is the value of 
the pion decay constant at zero quark mass. (Its normalization here is such that the 
physical value is /o — 93MeV.) In what follows we shall identify the quark mass niq in 
X with the PCAC quark mass rtiq^'^^ . According to the previous discussion this is a 
non-trivial choice because the lattice artifacts in (ratios of) the quark mass are rather 
different for arriq'-"^'-" then, for instance, for arriren- 

The sea quark mass dependence of the ratio of the pion decay constant in NLO of 
ChPT is: 

Rfss^^ = l + Ha--^)XR{NsLm + LR5)-^^aloga + 0{xji) ■ (2) 
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Here fss is the pion decay constant of a pion consisting of two sea quarks with mass 
Xs and fjiji is its value at some reference quark mass XR- is the number of mass- 
degenerate sea quarks (actually Ng = 2), L/j^. {k = 4,5,...) are Gasser-Leutwyler 
constants at the scale fi = fo^/XR the ratio of sea quark masses to the reference 
quark mass is 

a^^. (3) 

XR 

The analogous formula for the pion mass squares is: 



m 



2 

Rnss = ^^^ = l + 8{a-l)xR{2NsLm + 2Liis-NsLm-LR5) 

Note that instead of the scale dependent combinations (at Ng = 2) 

Lra5 = 2-Lij4 + Lij5 , -^^_R6845 = 4Lr6 + 2-Lij8 — 2Lr^ — Lr^ (5) 
one can also use the universal low energy scales A^^^ defined by jlUj 

A3 = 47r/o exp( -06845 ) , "6845 = 1287r^Lij6845 - 2 l°g Yevr^ 

A4 = 47r/oexp(a45/4) , 0^45 = 1287r Lji^^ + 2 log ^^^2 ' 

(6) 

The free parameters in Rfss and Rnss are XRi XR^Rib and XrLr6845- With 
the small number of points we have the linear fit with these parameters gives a good 
chi-square but relatively large errors: = 0.8 and 

XR = 30.8(9.4) , XrLra5 = 0.1398(86) , XRiR6845 = -0.0078(22) . (7) 

This corresponds to 

Lh45 = 4.5(1.1) • 10"3 , ^ = 23.3(8.2) , 

Jo 

LR6845 = -2.54(21) • 10-^ , ^ = 7.64(14) . (8) 

Jo 



Consistent results with smaller errors can be obtained if one takes the value of 
Xr. = 35.8(3.3) from the fit of the valence quark mass dependences (see next section) 
and performs two linear fits with the parameters xrLra^ and X-R-^-R6845) respectively. 
The resulting parameters are 

XrLra5 = 0.1443(15) , Lh45 = 4.03(37) • lO'^ , ^ = 21.4(1.5) , 

Jo 

XR^R6845 = -0.00896(86) , Lr6845 = -2.50(34) • 10-4 , ^ = 8.21(27) (9) 

Jo 
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Sea quark mass dependence of Rf 

1.2 
1.1 
1 

0.9 
0.8 
0.7 
0.6 
0.5 
0.4 
0.3 

0.2 0.4 0.6 0.8 1 

a 

Figure 1: Sea quark mass dependence of the pion decay constant. 
The straight dashed hne connects the first two points. 

and the fits are shown in figures Q and |2 

As these figures show, both Rfss and Rnss can be well fitted with the NLO ChPT 
formula. The fit parameters are within the expected range. For instance, the value 
of xr is rather close to the tree-level estimate -yH**"^"*^ / {r^f^f' ~ 40.3. (Here 

we used ro/o — 0.23.) The presence of a chiral logarithm which causes the curvature 
is clearly displayed in figure ^ where a straight line connecting the first two points is 
also shown. In Rnss the measured points are consistent with the presence of a chiral 
logarithm but the relative errors are large because all the values including the ChPT 
fit are very close to 1. This implies that the deviation from the tree-level behaviour 
m'gg oc xs is rather small. The results for the parameters in Q are close to the ones 
reported in jTj: the values for A4//0 practically coincide and the value of A3//0 is only 
slightly higher now. 

The extrapolated values of Rfss and Rnss at zero quark mass are, respectively: 
Rfo = 0.4228(60) , Rno = 1.0717(69) . (10) 

The value of Rfo togehter with Z^^af-jr from tableland tq = 0.5 fm imply for the 



16^x32 Xr = 35.8(3.3) 

|3 = 5.10 Lfi45 = 0.00403(37) 

K = 0.176 - 0.177 A4/fo= 21.4(1.5) 




- fit: 1 +4(a-1)xRLR45-XRaloga/167i^ 
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Sea quark mass dependence of Rn 
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Figure 2: Sea quark mass dependence of the pion mass-squared 
divided by the quark mass. 

pion decay constant at zero quark mass (/o) 

ZjVo = 121(5) MeV . (11) 

This result for Ng = 2 hght quarks compares weh with the phenomenological value 
/o = 93MeV if, as expected, Za = 0(1). 

3 Valence quark mass dependence 

We consider for fixed sea quark mass xs the valence quark mass dependence of /^r and 
as a function of the quark mass ratio 

e= — ■ (12) 
XS 

In our numerical data we determined the pseudoscalar mass and decay constant in rela- 
tively wide ranges of the valence quark mass ratios, typically ^ < ^ < 2. At the smaller 
quark masses {k = K2,3), however, for ^ < 1 exceptional gauge configurations appear 
which blow up the statistical errors and clearly influence the mean values themselves. 



16^x32 Xr = 35.8(3.3) 

P = 5.10 Lr6845 = -0.000250(34) 

" K = 0.176 -0.177 A3/fo = 8.21(27) 




fit: 1 + 8((j-1)XrLr6845 + XrO \oga/32n^ 
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Therefore, in most cases we restrict our fits to valence quark masses larger than the 
sea quark mass > 1). 

In the partially quenched situation several types of ratios can be constructed be- 
cause the pseudoscalar meson can be the bound state of two valence quarks (V^l^) and 
also a valence quark and a sea quark (VS). The PQChPT formulas for the ratios of 
decay constants are: 

Rfyy^^ = l+4(e-l)X5^55-^(l+01og^ 

+ DfvvxUC - l)+QfVVXsi^ - if + OixllogCxs) (13) 



and 

fss 



Rfys = ^ = i + 2{^-l)xsLs5 



+ lDfvvxUC-l)+QfVsxl{^-lf + 0{xllogC,Xs) ■ (14) 

The analogous formulas for the valence quark mass dependence of the (squared) pseu- 
doscalar meson masses are: 

2 

Rnvv = ^ = 1 + 8(e-l)x5(2L58-i>S5) + T-^[e-l + (2^-1) log 

^^^^ S S -LDiV^TT 



+ DnvvXsiC - 1) + QnVVXsiC " 1)' + 0{xs log t Xs) (15) 



and 



+ ^DnVVxU^ - 1) + QnVSxU^ " 1)' 

+ 0(x|logC,x|) • (16) 

In these formulas the Gasser-Leutwyler coefficients Lsk (A; = 4, 5, . . .) are defined at the 
scale fo^/xs and, in addition to the NLO terms, also the tree-graph (i.e. counterterm) 
contributions of the NNLO are included. Their general form is taken from ref. jllj and 
is discussed in more detail in section 2.1 of [J. The left-out terms of NNLO, which 
come from two-loop integrals, are generically denoted here by 0{x's^og(,). 

In addition to the single ratios Rfvv, Rfvs, Rnyv and Rnys it is useful to 
consider the so called double ratios which do not depend on any of the NLO coefficients 
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Lsk- The PQChPT formulas for the double ratios are: 

RRf^Jy^ = i + -^^(^_i_iog^) 



fvvfss 32iV,7r2 

+ Q/dx|(?-l)' + 0(x|loge,x|) (17) 



and 



+ OndxiK - 1)' + 0(,\slog«.Xs) ■ (18) 

In the PQChPT formulas (|13j) - (|18|) there are altogether 11 parameters. Three of 
them appear at NLO, namely with Ng = 2 

XR , XrLr5 , XrLr85 = XrC^Lrs - Lr5) (19) 

and the rest in NNLO: 

XRDfVV,nVV , XrQ fVV,fVS,fd,nVV,nVS,nd ■ (20) 

At the smallest quark mass fits with the NLO formulas are reasonable but for the larger 
quark masses the NNLO contributions are required unless the fits are restricted to a 
small range around = 1. 

An acceptable global fit with 11 parameters can be achieved if the valence quark 
mass dependence at all four sea quark masses is simultaneously considered. In this case 
one has to choose a reference sea quark mass XR take into account the relation 
between the NLO parameters 

Lsk = LRk - Ck log — (21) 
XR 

where the relevant constants are: 

C5 = Y^, C85=2C8-C5 = -^. (22) 

Fitting all six valence quark mass dependences {Rfw, Rfvs, RRf ■, Rnyv-, R'^vs-, 
RRn) there are reasonably good 11 parameter (linear) fits with — 200 ~ no. of 
degrees of freedom. A typical set of the resulting fit parameters is shown in table [21 

Comparing table El to the corresponding one (table 4) in ref. ^ one can see that 
most values are within statistical errors the same. This is also true for the NLO 
parameters defined at the scale 47r/o, which are now 

as = 1287r2LK5 + log = 2.06(42) , 
ass = 2as - as = 128tt''Lrs5 - log ^ = 0.583(45) . (23) 



Table 2: Values of best fit parameters for the valence quark mass 
dependence. Quantities directly used in the fitting procedure are in 
boldface. 



Xr 


35.8(3.3) 








0.1003(76) 




2.80(39) 


• 10" 


-3 




-0.0256(12) 


-TtOO 


-0.714(65) 


■ 10" 


-3 


■y'i.Dfw 

/\-Ji J V V 


-0.109(42) 


D fvv 

J V V 


-8.5(4.4) 


■ 10" 


-5 


V oQ f w 

^ix^ J ^ ^ 


-0.014(29) 


W fVV 


-1.1(2.3) 


■ 10" 


-5 


x'rrQ fVS 


-0.0177(94) 


(•J fVS 


-1.39(81) 


■ 10" 


-5 


x%Qfd 


-0.0180(31) 


(J fti 


-1.41(13) 


■ 10" 


-5 


XjlDnVV 


-0.134(21) 


DnVV 


-10.46(93) 


• 10" 


-5 


X%QnVV 


-0.087(13) 


QnVV 


-6.77(30) 


■ 10" 


-5 


X%QnVS 


-0.0394(44) 


QnVS 


-3.07(24) 


■ 10" 


-5 


XjlQnd 


0.0077(48) 


Qnd 


0.60(26) 


■ 10" 


-5 



The value of as is practically the same as in table 5 of ^ whereas ags is slightly smaller 
now. 

The tree-graph NNLO contributions play an important role in the global fits of the 
valence quark dependences, especially at the two larger sea quark mass {n = kq and 
K = Ki). At the two smaller sea quark masses NNLO is substantially less important. 
This is illustrated in figure |3 where the 11 parameter fit for Rfvv is shown together 
with the NLO contributions alone. 



3.1 0{a) terms 

The fits above have been performed with the continuum formulas - without 0{a) or 
any other lattice artifacts. The fits are reasonably good and the resulting parameters 
are quite similar to those obtained in ref. ^ where the 0{a) terms have been taken 
into account in the (PQ)ChPT Lagrangian according to ref. J2|- It has been observed 
already in [1] that the parameter in the chiral Lagrangian characterizing the magnitude 
of 0{a) effects 

_ 2Woacsw 

P = 72 (24) 

JO 

is rather small compared to the quark mass parameter x in ©• Fitting the ratio 
rjs = ps/xs separately for the individual sea quark mass values we obtained increasing 
values for increasing sea quark masses: 0.02 < r]s < 0.07. 

The parameter p should be independent of the quark mass because the quark masses 
are the other expansion parameters in the chiral Lagrangian. This means that a quark 
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Valence quark mass dependence of RfVV Valence quark mass dependence of RfVV 




Figure 3: The linear fit for Rfvv for different sea quark masses. 
The NLO contributions alone are shown by dashed lines. 



mass dependent ps incorporates some higher order effects proportional to some power of 
aruq. (For instance, a linearly incrasing value of 775 corresponds to ps oc {amqf' .) Since 
the observed values of p are small anyway it is interesting to consider the behaviour of 
the chi-square as a function of p if the linear fits are performed for fixed p . Because of 
the presence of another new parameter describing 0{a) effects in the chiral Lagrangian, 
the linear fit has 12 parameters for p ^ (instead of 11 for p = 0). As it is shown by 
figure 13 the of the fit has a minimum near p = r] = and becomes extremely large 
already at |r/| ~ 0.1 where the absolute value of p is 10% of the value of the reference 
quark mass parameter XR- 

Another way to investigate the importance of 0{a) effects in our data is to consider 
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p dependence of chi-square of a linear fit 
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Figure 4: of the linear fit of the ratios Rfvv, Rfvs, RRf, 
Rnvv, Rnvs and RRn for all four sea quark masses at fixed values 
of the 0{a) parameter p in the chiral Lagrangian. 



the following combination of double ratios: 

RRn + 2RRf-3 = {Qnd + 2Qfd)xUC-lf + 0{xllogC,Xs) • (25) 

As this formula shows, this combination vanishes in next-to-leading order and only 
NNLO and higher orders contribute to it. On the lattice there could also be 0{a) 
contributions which can be parametrized as 

(£_1)2 (£-1)2 

RRn + 2 RRf -3 = 16pLsAW6-^ — TT " P ' 



[2(i-e2) + ioge + 3£2ioge] (e-i-^iogo 



327r2£(e + 1) 



327r2£ 



+ 0{p'',x^' 



(26) 



Here only the linear piece of the r]s = /o/x5-dependence is kept because r]s is small. 
-^541^6 = — Ws6 is a new parameter appearing in the 0{a) terms of the chiral 
Lagrangian |12| . 
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The linear fits with XsQn2f = XsiQnd + '^Qfd) in and with p in (|26|) . respec- 
tively, are shown in case of the smallest sea quark mass {k = K3) by figure |SJ As 
this figure shows, the NNLO fit with x'sQn2f is better (x^ = 1.3) than the one with 
the leading 0{a) term proportional to p (x^ = 7.2). For simplicity, the parameters 
Xs = 11-7 and Ls^we = 0.001 are fixed in this latter case but taking other values does 
not change the qualitative picture. 

Valence quark mass dependence of RRn+2RRf 
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w -0.002 
c 

tr 

^ -0.004 
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-0.008 
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Figure 5: Comparing the NNLO fit (full line) with the leading 
0(a) fit (dashed line) for {RRn + 2RRf - 3) at k = K3. 

At the larger sea quark mass values the fits with the leading 0{a) terms behave 
similarly to figure |SJ This supports the fact that the 0{a) terms are not important in 
our numerical data. As shown by figure 31 good fits can only be obtained at rather 
small values of 77 = p/x- In contrast, the NNLO contributions are very important 
especially at our larger sea quark masses. 

4 Studies of the updating algorithm 

The numerical simulations have been performed by the two-step multi-boson (TSMB) 
algorithm This dynamical fermion update algorithm is based on the multi-boson 



1 

I 

T 


1 1 1 1 1 
1 6^ X 32 lattice 
P = 5.10 
K = 0.177 

Xs^Qn2f= -0.0068(29) 
Tls=p/Xs = 0.041 (21) 


1 


1 1 1 1 1 
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representation of the fermion determinant 14 and in its present form it incorporates 
several modern ideas of fermionic updating: the global correction step in the update 
jl5j . the final reweighting correction ^HI and the determinant breakup II 7j . 

Our error analysis is based on measuring the autocorrelations of the quatities in 
question ^l^j therefore we can estimate the computation cost based on the integrated 
autocorrelations Tint- In our previous papers 012011^ we proposed an approximate 
formula for the cost 

a„,c^F(am,)"2j7 (27) 

where arug is the quark mass in lattice units and the number of lattice points. The 
overall factor F depends on the quantity under investigation. If we count the cost 
in terms of the number of floating point operations necessary to perform an update 
sequence with length Tint then the present simulations on 16^ • 32 lattice are consistent 
with 

Fplaguette ^ 7 • 10^ , ~ 10^ , Ff^ < 4 ■ 10^ . (28) 

In case of f-j^ we only have an upper limit on Tint because the gauge configurations stored 
for the measurements were statistically independent. These numbers are somewhat 
smaller than our previous estimates in [201 which is due to a better tuning of 
algorithmic parameters. In particular, these simulations were done with a determinant 
breakup A'^^ = 4 which means that the fermion determinant of the two degenerate 
flavours {Nf = 2) were reproduced by 4® (A'^^ = ^) flavours. Another important point 
is the frequent call of the global heatbath update of the multi-boson flelds which every 
time gives a statistically independent boson conflguration. 

If we take the plaquette expectation value as the worst case, then at the present 
quark masses and lattice spacing this cost estimate is similar to previous estimates 
(see, for instance, the formula of A. Ukawa j22]) but considering the more interesting 
cases of m-,^ or 7,^ there is a substantial improvement by an order of magnitude or 
more. In addition, towards large volumes, smaller quark masses and/or smaller lattice 
spacings the scaling of the cost estimate in is better: for flxed lattice spacing the 
cost increases as m~^il and decreasing the lattice spacing and keeping the physical 
parameters flxed the cost behaves as a^^. This has to be compared to the estimated 
behaviour in j22j m~^il^/^ and a"'', respectively. 

5 Summary 

The quark mass dependence of the pseudoscalar mass and decay constant in our nu- 
merical data can be well fltted with the continuum (PQ)ChPT formulas. It has been 
already observed on the 16^^ lattice in ref. ^ that the 0{a) lattice artifacts at our 
gauge coupling (3 = 5.1, corresponding to a lattice spacing a ~ 0.2 fm, are small and 
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one can obtain reasonable fits by omitting them. This conclusion is strengthened by 
the new 16'^ -32 data and therefore here we based our estimates of the chiral Lagrangian 
parameters on fits with the continuum formulas. 

The use of the ratios of the PCAC quark mass as the variable in comparing the 
simulation data to chiral perturbation theory is essential. Taking other quark mass 
definitions, for instance Hren = ZmifJ-o — fJ-cr), would be the source of large lattice 
artifacts at our lattice spacing. 

The sea quark mass dependence of fj^ and /mg can be well described in our quark 
mass range 0.855 < Mr < 2.13, which roughly corresponds to \ms < mq < |ms, by 
the NLO ChPT formulas. The obtained estimates of the relevant Gasser-Leutwyler 
constants are, according to section |2| 

^ = 8.21(27) , ^ = 21.4(1.5) . (29) 

Jo Jo 

The functional dependence of the ratio of as a function of the ratio of quark masses 
clearly shows the presence of chiral logarithms (see figure ^ . This observation is in 
agreement with the results in a recent paper of the UKQCD Collaboration which 
came out during writing up this paper. 

In the valence quark mass dependence of the same quantities, in addition to the 
NLO terms, the higher order NNLO contributions appear to be important - especially 
at our two larger sea quark masses. But, as shown by figure |1J the importance of the 
NNLO terms is considerably reduced at the two lighter sea quark masses. Our best 
estimates for the relevant Gasser-Leutwyler constants at the scale 47r/o are, according 
to (EHJ: 

05 = 2.06(42) , 2a8 - as = 0.583(45) . (30) 

The errors quoted in (I^U]) and (|5n|) are only the statistical ones. In order to decrease 
the systematic errors simulations at still smaller sea quark masses would be useful. 
Since our lattice volume is relatively large (L ~ 3fm), finite volume effects can be 
expected to be small (see |24| I25j). For the moment we have no direct handle on 
the magnitude of the remaining non-zero lattice spacing effects. These should be 
determined by performing simulations at smaller lattice spacings. 
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